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ABSTRACT 


In this report a technique has been proposed to accurately estimate the 
Doppler frequency and demodulate the digitally encoded telemetry signal that 
contains the measurements from the balloon instruments. Since the data are 
prerecorded, one can take advantage of noncausal estimators that are both simpler 
and more computationally efficient than the usual closed-loop or real-time esti- 
mators for signal detection and carrier tracking. Algorithms for carrier fre- 
quency estimation, subcarrier demodulation, bit and frame syrchroniza*- : < • are 
described. A Viterbi decoder algorithm using a branch indexing technique has 
been devised to decode constraint length b, rate 1/2 convolutional code that is 
being used by the balloon transmitter. These algorithms are memory efficient and 
can be implemented on microcomputer systems. 
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I. 


INTRODUCTION 


The Venus Balloon (Vega) Mission Is a joint French-Soviet-and US project to 
place two balloons on the atmosphere of Venus. The balloons will be tracked by 
Very Long Baseline Interferometry (VI.BI) to determine their positions and veloci- 
ties on the disc of Venus. Besides this, the gondola contains an instrumentation 
package that sends telemetry data directly to earth using an L-band transmitter. 
This mission is described more fully in reference [1]. 

The scientific data recovery of the Venus Balloon mission requires accurate 

estimates of the Doppler frequency of the balloon transmission as well as the 
demodulation of the digitally encoded telemetry signal that contains the measure- 
ments from the balloon instruments. Toward this end, the balloon signal has been 
designed to provide an easy method to determine adequate starting parameters from 
this weak signal. Specifically, an unmodulated carrier wave is transmitted at 
the beginning and end of the data transmissions for a period of 30 seconds. 

A known bit sequence follows the carrier wave such that the starting time 
can be estimated accurately. To aid in carrier tracking through the frame, 

43° modulation is used such that half the transmitted power remains in the 
carrier wave. The subcarrier wave is separated far enough from the carrier that 
the information spectrum does not approach the carrier frequency; thus, it is 
easy to filter the carrier out of the signal for accurate Doppler extraction. 

The Doppler shift of the balloon transmission depends upon the balloon's 
location in the atmosphere and upon the motion of the balloon within the atmos- 
phere. As the balloon drifts across the disc of Venus, the Doppler shift should 
drift from high to low frequency. The total drift should be constrained within 
2 kHz. The width of the modulation spectrum requires an extra 500 Hz. Thus a 
total spectral width of 2.5 kHz needs to be preserved if there is no a priori 
knowledge of the initial frequency. As there is some a priori knowledge, probably 
it is safe to preserve 2 kHz which requires that roughly 4k samples per second 
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m n 


be recorded. Eight or twelve bit samples are fully adequate as no appreciable 

quantizing errors will exist, since the signal to noise ratio in the 2 kHz 

bandwidth is well below unity. Telemetry transmissions occur at a minimum period 

9 

of one half hour. The total recorded telemetry data represents less than 10 
unprocessed bits. 

The signal to noise ratio is given by 


P T C , 

= ~ T nil y -- /kTB 
4ttD" 


where 

P 


min 


'R 


G 

S 

A i 

D 

k 
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transmitter power in watts 
minimum transmitting antenna gain 

the effective receiving antenna collecting area in sq. meter 
the distance in meters 
Boltzmann's constant in watts/K/Hz 
the effective receiver system temperature in K 
the receiver bandwidth in Hz 
where the full 2 kHz bandwidth is used. After filtering to a 1 Hz bandwidth, 
the signal co noise ratio increases to roughly 63. In this narrow bandwidth, 
the signal should be clearly visible on a spect r um by spectrum basis. Tor 
example, using DSS-14 with an l.-band FET amplifier would give P /P = 3.16 x 
10 ^ where 

P T - 5 
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1.63 x 10 3 
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1.38 x 10 
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II. INITIAL DETECTION TECHNIQUE 

a. Determination of Filter Passband Shape 

Prior to the arrival of the actual Balloon signal and periodically while no 
signal exists, the receiver passband 3hould be measured accurately. Since the 
passband shape is determined by the offset baseband filter and any aliasing due 
to the discrete sampling, it is not expected that the shape will vary signifi- 
cantly over any time scale of interest to the detection of a single telemetry 
frame. An accurate enough passband spectrum can be obtained by averaging 200 sig- 
nal free power spectra to get roughly 7% statistics. If each spectrum uses 4096 
samples, the time to acquire this average noise spectrum Is 200 x ^ - q - q - q = 204.8 s or 
about 3.4 minutes. 

Procedure : 

Using the signal free region of the recorded data, perform the Real Fourier 

Transform (RFT) algorithm. (Details of the theory of the RFT have been supplied 

as NASA Tech Brief No. NPO-11649.) The real Fourier transform results in place 

of the data array that contained 4096 points. Next form the power spectrum from 

these data. The power spectrum will result in N/2 + 1 points (2049) as follows: 

The first point in the transformed array is the a term, the second point is the 

a , n term, the remaining N-2 points are complex, thus a = d + j d , 

N/2 1 J ** 

*2 • d 5 + i V "k “ d 2k+l + j <*2k« •••• a N/2-l ’ Vl + •> V a 0 * d l + J °’ 

a N/2 = ^2 + 3 ^ w ' ierc u k a complex component of the transformed array and d^ 
are the result of the RFT program. 

To begin, clear the power array which contains N / 2-4-1 points. Next 
accumulate 200 spectra points to this array by forming the power from the array, 
where 

P lk) «- P (k) + a. a* ; k = 0 N/2 

n n k k 
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Save Che P array to use as a baseline reference to detect the unmodulated 
n 

carrier. Note that the standard deviation of the noise spectrum will be about 
77 . of the mean value. 

b. Baseline Removal 

Spectra measured during the carrier portion should contain a strong line at 
the carrier frequency. Due to acceleration on the balloon, the carrier frequency 
could vary by as much as 0.5 Hz/s but rather uniformly for periods as long as the 
data transmission period. Thus it is not advisable to accumulate the spectra for 
any appreciable period as smearing would result. The signal power spectra should 
be determined in the same way as the noise power spectra. Let P g be a signal 
spectrum resulting from a single RFT. Then the effects of the passband shape 
can be removed by division. 

P (k) «- P (k)/P (k) k - 0 N/2 

s s n 

c . Threshold Detection 

The object of this section is to provide a way to automatically search for 
the spectral line and to determine if a spectral line has been detected in the 
present data frame. Let k^ and k^ represent the limits of the search where k^ 
could start at zero and k could range to N/2. Let there be M spectra containing 

M 

pure carrier. In reality a single spectrum is formed from 1.02 s of data, and 
there are 30 s of pure carrier, so there will be roughly 29 such spectra. Let 
L(n) , 1 < n £ M be an array that stores the location of the detected 
line for each spectrum. 
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<cedure: 


Since each spectrum contains noise having statistics of Chi , and some 

point will always be maximum, a threshold must be established that 

assures that the signal power is significantly larger than the maximum 

noise power. This can be done only within some confidence level since 

the noise statistics are not limited in range. Note that the signal is 

normally divided between two spectral points. Search the normalized power 

spectrum for the max power. This point is then P (k ) ■ P . Now form 
r r r s m max 

the following statistics for each spectrum: 

N/2 

nTT Z P . <k) 

k*0 » k^k 

m 

N/2 

-S-/2- Z ‘ >2 S <k) ' ^ 

k-0 , k*k 

m 



Test P > M + 9a 
max 

If the above is true, this is a tentative detection. If not true 

continue to the next spectrum. To be statistically certain of a 

detection, save the value of k and process the next two spectra. If 

m 

all three yield a value of k within a range of two, the detection is 

in 

reasonably assured. This can be checked by finding the maximum and 

minimum values of k . If the difference between the max and min is 1 

m 

or 0, the detection is established. Record the starting sample number 
associated with the first point of the first spectrum and maintain it 
as a starting reference point for further processing. 


5 



2. Beginning at this time reference, use the above procedute to fill the 

array L(n) with the successive values of k • Process 29 spectra. 

m 

d. Carrier Frequency Estimation 

The array L(n) contains the trace of the cw line as a function of time. Tt 

is expected that this will form a straight line having a small slope. There is 

some probability that bad estimates of k exist. It is likely that these will be 

in 

widely separated from the main sequence of points. Thus it is necessary to 
attempt to eliminate the bad points if they exist. 

Procedure : 

1. Find the mean of L. and mark all points deviating from the mean by 
more than 5 by replacinp the value of L(n) with -1. 

2. Find a line of regression by linear least squares that represents L(n). 
Skip all points containing -1 in this analysis. The result of this 
analysis is an equation that gives L ■ a x + b. In the above analysis, 
the values of x should be the sample number. For example, x^ should 
be associated with the center sample number for the first spectrum. 

That is, x^ is sample number 2048 counting from the sample reference 
number. In this way the variation of the instantaneous frequency can 
be predicted relative to the reference sample number. Thus if i is 
the sample number and i ^ is the reference sample number, the fre- 
quency (in frequency number) is given by 

L ■ a (i - i J+b 
ref 

and 

X “ (1 ref • 1} + f + n N 
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where n is the spectrum number and L is a measure of frequency cell 
number as a function of sample number. It is useful to know both 
phase and frequency as a function of sample number where real time is 
referenced to the first sample (i ^). If S R is the sample rate and 
N the number of samples in the RFT, then 

S R 

Let w(i) = 2tt — [(i - i ref ) a + b] 


t = (i - 


W /S R 


Then 


w ( t ) = [a S R t + b] 

Sr 2 

t(t) = 2 "y U S R t /2 + bt] 


Thus 


w = 


S R 


and w = 


tS r 3 

N 


3. Repeat the above procedure for the end-of-frame cw spectra. This 

forms a second equation of frequency versus sample number. Let be 
the first and L the second, then some measure of frequency stability 

D 

over the data frame can be determined by evaluating the difference of 

L. and L_ at a point near the center of the frame. If this difference 
A B 

is too large, it may be difficult to track the carrier through the 
frame . 
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III. DOPPLER RECOVERY AND THE EQUATION OF TIME 

a. Filtering and Data Reduction 

All balloon motion relative to the center of Venus is unmodeled in the 
frequency drifting of the first local oscillator. Thus it is important to know 
the Doppler shift accurately to know the radial velocity and acceleration of the 
balloon as well as to be able to demodulate the carrier wave. L A and Lg give a 
first order prediction of that Doppler, but furtner prediction may be carried 
out by filtering out the sidebands and tracking the carrier through the frame. 

The frame consists of 330 seconds of data, and there are 1.32 x 10^ samples 
through the frame. Since this is far too many to store in the computer memory, 
it is important to perform some data reduction. We are interested in frequency 
departures relative to our initial frequency line of regression. Normally, 
departures of more than a few Hz from this line are unexpected (note each 
frequency bin has a width of 4000/4096 = 0.9765 Hz). To get some bound upon 
this, determine the maximum difference between |L - L j over the whole time 

A D 

span. The width of the filter can be set at twice this value. Call this 
bandwidth B. Note that the Nyquist sampling requires only 2B samples per second. 
Since B will be only a few Hz, only four or five samples per second are required 
to specify the carrier signal. Thus, roughly two thousand memory locations 
are needed to store the information for the entire frame. This data reduction is 
achieved by heterodyning the input signal ^samples) with a reference carrier 
based on L A and filtering. The filter can be a simple running average of 1/2B 
seconds of data; however, it is convenient to form complex samples at a rate of 
1/B using a complex mixer. Figure 1 below shows the block diagram for a complex 
mixer. 
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Figure 1. Baseband Filter for Carrier Separation 


Procedure : 

1. Determine B as discussed above. 

2. Let T = 1/B. The number of samples to average in the filter is then 
4000 T = M. 

3. Beginning at the reference sample form the products 

w 2 

S^^ Cos [2ir(wt + | t )] t is time 

and • 2 

Sin[2ir(wt + y t ) ] 

Then produce the average of M such samples such that 


< k+1 >! + ‘ref 


C, - 


i=(k-l)^r + i 


Si Cos [2r(wt + 
M 


t 2 > 


ref 


~ J 


S. Sin [2 


w 2 

(wt + f t ) ] 


4. 


Save C, 
k 


for further analysis. 
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b. Estimation of Doppler Model 

In this section, we attempt to construct an accurate model of the carrier 
wave to be used as a phase coherent reference for demodulation of the subcarrier 
signal. This is essential to prevent the true carrier and the estimated local 
carrier from forming a beat that modulates the subcarrier amplitude. The phase 
coherent local oscillator is constructed from the carrier samples C^. The extent 
to which the carrier model can represent the data depends upon many factors 
such as the signal to noise ratio, the phase stability of the transmitter, the 
motion of the balloon, the atmospheric turbulence and many other stochastic 
errors and unmodeled effects. Thus it is not possible to know a priori the level 
co which coherence can be established. For that reason, it is necessary to know 
when the carrier tracking loop is locked as tightly as possible. A simple 
criterion is to establish whether the estimated carrier signal has fallen signi- 
ficantly below the estimated incoherent power level. Thus we attempt to lock the 
carrier model to the carrier samples as tightly as possible by narrowing the 
effective bandwidth until the signal level begins to drop. The effective band- 
width is controlled by the duration of the data span used in the processing. 
Increasing the span to 30s should be possible most of the time unless the balloon 
is undergoing severe turbulence or accelerations. 

Procedure : 

Nonlinear Least Squares Regression 

Let the carrier signal model be given by 

j ($ + wt k - w/2t£) 

\ =Ae 

where 

A is an amplitude coefficient 
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^ is an initial phase reference 

t^ is the time of the kth sample and k is centered in the span of data 
w is the frequency 

w is the rate of change of frequency 

Let the error between the measured samples and the model be given by 

N 

c ■ 2 H k l c k - “kl 2 

k*-N 


where is a set of weights that are normally set equal to the reciprocal of the 
noise variance, and 2N+1 is the number of samples in 30 seconds. In this case 
the weights can be set equal to unity unless the true covariance matrix is needed. 
We desire to minimize e by adjusting A, <J>, w, and w. Let x represent the vector 
A, <p, w, w and Ax represent a vector of small corrections to x. 

Ax * [AA, <p. Aw Aw] 


Then 


9M R j(<|> + wt R + w/2t£) 

~ = 3x^ = 6 


d \ 


9M. 

9x 


k j(4> + wt fe + w/2t k ) 

= j A e 



3 \ 

9x 


j(4> + wt + 


TT-JAt. 


“ /2t k> 


3\ 3M , J(* + wt + w/2t 2 ) 

■ 5T ’ J A c k /2e 
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A minimum In e requires that 



1 


1, 


,4 


Thus 


- L 2R eK <( > -V S' 0 ‘-I.-.* 

k— N 1 


may be expanded about the initial value of x. If we substitute only 
the linear terms of the series for M^, we arrive at the normal equations. Let 




M + 

ok 2mJ 5x . 

j-I 


Ax , 


! X “X . 

I j °j 


Then 


N 


2 R 


k— N 


* 

4 N 1 


3M k 

W, (C. -M . ) *7"“ 
k k ok 3x^ 

| j-1 k— N 


w 

k 3Xj 3x^ 


Ax, 


3M k 

evaluate -r-^ at x ■ x 

3x j 


The normal equation forms four equations with four unknowns where Ax^ are 
the unknowns. This solution must be Iterated until the Ax^ are suitably small. 
Thus the next estimate of x q ■ x q + Ax. The meaning of suitably small is either 
defined from the limits of the measurements as defined by the covariance matrix 
(the inverse of the normal equation matrix) or by what is reasonable. Reasonable 
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in this case is determined by theoretical limits or by how much error we can stand. 
In general a phase error of 10° is about as large as acceptable. This constrains 
A$ < 0.17, Aw < 0.01 and Aw < 0.0015 for a 30s span. In general, the iteration 
should be driven to about a tenth of these values. 

Accurate starting estimates are usually required to begin the nonlinear least 
squares regression analysis. The prediction equation given by L A should normally 
determine the frequency to an accuracy slightly better than 0.5 Hz. In such a 
case the maximum integrated phase error approaches 360® in two seconds. The least 
squares procedure requires that the phase error be less than 180° and preferably 
less than 90° over the data span if convergence is to be assured. Assuring the 
phase error is small enough either requires that the preliminary estimates of w 
and w are adequate for the time span considered qr that the time span is short 
enough. If the time span is too small, the signal to noise ratio will be too 
small for adequate detection. Since most of the w term has been heterodyned out 
by the model, it normally can be initialized to zero. Determination of w can 
be done in several ways, each of which has merit under certain conditions. We 
suggest the following possibilities: 

1. Under normal conditions, the NLLSR can be used with a short span of 
data (about Is). Once convergence is established, the span can be 
doubled using the new estimate of w. This can be continued until the 
30s span either converges or the estimate of A begins to drop, indicating 
inadequacy in the model or phase incoherence of the signal. 

2. A sequence of Is estimates can be made over the 30s data span using 
NLLSR. The resulting w then can be used to form a linear equation of 
frequency giving w and w as starting parameters. This procedure has 
been tested and works well if the signal to noise ratio is adequate. 
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3. This last procedure may be useful if convergence cannot be obtained 
using 1 and 2 above. Use the complex fast Fourier transform (FFT) to 
form the power spectrum of the C^'s. The span of data could extend to 
30s giving 1/30 Hz resolution. A peak detection scheme similar to that 
of Section lie can be used to find the max-timim power and its location. 
This procedure will not work if w is too large* but does work for 
smaller signal to noise ratios than procedure 2. 

A first estimate of Ae^ can be formed by simple correlation of the model 

with the data C . 

k 


A 

o 


1 

2N+1 


N 


£ 

k*-N 


-J (ut k + 2 C k> 
e 


Thus 



a ^ ^ 

<p =* Arctan (Im A/Re A) 

As we progress through the data, it is best to determine the starting values from 
the results of the previous data. However, A and 4 can always be computed quickly 
from the procedure above. 

Experimentation with some carrier models has shown that some procedure for 

~2 

detecting nonconvergence is necessary. A suitable estimate is to assure that A 
does not fall significantly below the estimated carrier power level. An estimate 
of the incoherent power is given by P where 

N 2 

* = 2N+I 2 ^ C kl " * 

k*-N 

N is an estimate of the noise power which can be formed from the original signal 
statistics by noting that the power in the 2 kHz band is mostly noise. Thus the 
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value of N can be related simply to the statistics of by noting that M samples 
of the original signal were added together and the gain of the complex mixer is 
1/2 per side. 

K 

*21 V ^ 9 

N ■ « "s • M 2isr 2 - S 1 

i--K 

2 

Here a is the variance of the samples, S. in the data span currently being 

S 1 

analyzed. This assumes are zero mean and nearly white. 

The results of the above analysis provide a series of equations for estimat- 
ing <j>(t) that are valid over a 30 second strip of data. The 30 second spans are 
formed every 15 seconds. These data are intrinsically valuable as scientific 
output and are used for the carrier demodulation discussed in the next section. 
Clearly, the previous operations provide a great data reduction in specifying 
the phase function since each 30 second strip of data is represented by four 
numbers and there are only 22 such strips over the data frame. 

The above procedure provides no guarantee that the phase is continuous across 
the thirty second boundaries. However, since the equations are valid over a 15 
second period overlapping consecutive strips, it is possible to force the phase to 
be continuous by using a weighted combination of two adjacent equations to form 
the phase function. Figure 2 below shows the form of the weights. 
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Figure 2. Linear Weights for Interpolation of Phase 


The Interpolation Is formed as shown In Figure 2 where 


W k (t x* At ’ " fc k+l ” t k 


and 


where 


♦ (t x ) - ( 1 - w k Ct x ))^ k (t x ) + VSc^k+l^x 5 


w k 2 

$. (t ) - *. t + t 
Y k x k x 2 x 


and w^w^ are from least squares regression. Thus the continuity of phase across 
the boundaries Is assured. 
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IV. CARRIER DEMODULATION 


The carrier demodulation ia eaaily carried out uaing the phaae function of 
the previoua aection. If we presume that the binary (square wave modulation) ia 
B(t) , then the modulated signal is of the form 

J[rr/4 B(t) + 8(t)J 

S(t) ■ S e 

where B(t) ■ +1 and S is a real amplitude weight. 8(t) is presumed to be composed 

A 

of drifts in the transmitter and Doppler shift due to balloon motion. 9(t) esti- 
mates both phase components , thus the recovered modulation r(t) Is given by 

-j 8(t) 

r(t) - S(t) e 

^ ^ y 2 

e(t) • ♦(t) + w A t + y t 

♦Note 1 


When 8. - 0. Is small as it should be, 
Q 0 


r(t) * S ~ + j sin | B(t) 


Note that verification of proper phase demodulation can be checked easily by 
simply averaging r(t). Since B(t) is zero mean, only the real part of r(t) 

should remain finite. Note that the modulation is contained in the Imaginary 

/2 

part of r(t) and has an amplitude of — ^ s - 
♦Note 1 

S(t) is the original sample set beginning 12 seconds after initial detection, 
thus this is 48000 samples past i^. 8(t) is then defined at the total phase 
at that time. t(t) must be composed of the original phase due to and the 
departures from given by $(t) where w^ and w^ are from (See page 7). 
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V. SUBCARRIER DEMODULATION 


First note that between times 31-3/8 sec and 34-1/4 sec we have an all zero 
sequence. Therefore there are no transitions in symbols. (Look at the output 
symbol >f quence of convolutional code for frame sync inputs.) Therefore we can 
demodulate the subcarrier and estimate the phase of subcarrier as shown in Fig. 3. 
In Figure 4 the output samples of the integrators in subcarrier demodulator are 
shown as a function of T 

sc 


Sin ( W|c tj) 



Fig. 3 Subcarrier Demodulator 




Fig. 4 Output Samples of Integrators in Subcarrier Demodulator vs. T 
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IF B 0 and C _< 0 
IF B £ 0 and C 0 
IF B < 0 and C > 0 


T x C 
SC 

T SC “ 4 (C-B) 

T sc (c ♦ 2B) 

T SC “ 4(C + B) 

T SC (3C - 2B) 
T SC “ 4(C - b") 


IF B > 0 and C > 0 T<!r - T SC (3C * 4B) 

SL A (C + B) 

after finding t we can shift the whole subcarrier signal by Tg C (delay by T^ c >. 
Then we can demodulate the. whole data frame by the subcarrier to get data samples 


dU^. 


VI. BIT SYNCHRONIZATION 

The bit synchronizer is shown in Fig. 5. 



Fig. 5 Bit Synchronizer 


Let d ( t ) be data samples. Suppose during one symbol time there are N samples 
d(tj). In upper branch of Fig. 5, we sum the N samples d(t^) starting at time 
kN + 1 + ‘ (- ^ j) . Next we detect the sign of symbol by passing through 
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a sign (*) function 


{ 1 ; x > 0 
-1 ; x < 0 


Next wa dAttet thA transition In consecutive symbols a k by finding onA-half of 
dlffarancA between consacutlva datactad symbol signs. If - 0 means thara Is 
no transitions, I k ■ + 1 naans thara Is a transition. If I k ■ + 1 then wa hava 
symbol transition fron + 1 to - 1 and if 1^ - - 1, than we have symbol transition 
from - 1 to ♦ 1. 


u 

In lower branch in Fig. 5 wa sum N samples d(t 1 ) starting at time kN - y ♦ 1 
+ 6. The result is sample Y k . Ua multiply Y R by I k to get Z R . We repeat this 
operation M times (M * 60), starting at k - k Q where time k Q N is roughly at the 
beginning of (30 - 42 sec) interval. Next wa find the result of the average 
of all M samples as c. Mow if c is positive wa increase 6; if c is negative 
we decrease 5. We repast the whole operation until we get smallest c possible. 


VII. FRAME SYNCHRONIZATION 


We get the smallest possible c from bit synch. We store all detected 
symbols a k in a shift register. Next we compare the pattern of symbols that we 
have at the output of encoder with contents of shift register. We shift this 
pattern until we match the pattern with the content of shift register with 
minimum number of discrepancies. Minimizing 0 in Fig. 6 can result in frame 

A A 

sync, where b k c(0,l) is related to a k by 


b 


k 



a 


k 



is exclusive OR. 
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Fig. 6 Frame Synchronizer 


The input frame synchronization bits to the convolutional code during the 
time interval (30-42) seconds are 

17(0), 5(1), 3(0), 2(1), 1(0), 3(1), 1(0), 1(1), 1(0), 1(1), 4(0), 1(1), 

2 ( 0 ), 1 ( 1 ), 1 ( 0 ), 2 ( 1 ), 2 ( 0 ). 

The corresponding output pattern of frame synchronization symbols of convolutional 
code is 

10 (X) , 24(0), 2(1), 2(0), 1(1), 1(0), 1(1), 5(0), 2(1), 1(0), 2(1), 1(0) 

2(1), 1(0), 1(1), 1(0), 2(1), 1(0), 3(1), 1(0), 3(1), 7(0), 3(1), 

3(0), 3(1), 1(0), 2(1), 1(0), 2(1), 3(0), 2(1), 1(0), 1(1). 

Where the number in front of parentheses indicates the number of consecutive 

bits shown in parentheses, the X means unknown bit. The contents of encoder 
at time 42 is: 



INITIAL STATK 
Of tNCODKR 
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VIII. VITERBI DECODER 


This section explains the Viterbi decoder algorithm for a constraint length 
K * 6, code rate r * 1/2 convolutional code. The code is shown to be a trans- 
parent convolutional COde.^ The weight distribution of the code has been computed 
by numerical analysis. The code state table, received sample quantization, 
metric quantization, and Viterbi algorithm using a branch indexing technique is 
explained in detail. The theoretical bit error performance of the code and the 
simulation results of this Viterbi algorithm is given. 

a. Code Structure 

Con? ’er a K = 6, r = 1/2 convolutional encoder which is illustrated in 
Figure 7. 



Figure 7. Code Structure 


^The code has been provided by mission as in Fig. 7 without any further informa- 
tion. The description of code was not available in the literature. 
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This code is a transparent code, i.e., the l's complement of the input bit 
sequence results in the l's complement of the output bits. To show this, note 
that the output code symbols can be represented in terms of the input bits as 


X lk - U k © U k-1 © U k-2 © U k-4 © U k-5 
^k * U k © U k-1 © U k-5 

where represents the "Exclusive OR" operation. Let IL denote the l's com- 

i » 

plenent of IL . Now if we replace by IL and Xy by X_.j ; i=l,2 ; j = k, 
k-l,....k-5 in the above code input-output relations and if the equalities hold, 
then the code is transparent. This can be shown by using the following relations 


U j * U.@l j = k,k-l,..k-5 

x y * ’'ij © 1 1 * 1 - 2 

in the above code input-output relations. When the code is transparent we can 
use outer differential encoding and decoding for inner channel coding t j resolve 
the 180° phase ambiguity. 

The code weight distribution has been found by first finding the transfer 
function bound [2] for the code and then using a technique in f 3 ] to get the 
weight distribution. The minimum distance of code is 8. The code weight distri- 
bution is shown in Table 1. 

This code differs from the K = 6, r = 1/2 Odenwalder code [4], since it is 
transparent and the structure is different. The performance is slightly in- 
ferior to the Odenwalder code. 
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Table 1. Code Weight Distribution 


Weight 

d 

Number of 
Adversaries 

Number of Bit Errors 
in the Adversaries 

a d 

8 

3 

6 

9 

0 

0 

10 

13 

60 

11 

0 

0 

12 

~ 80 

469 

13 

0 

0 

14 

- 

3340 


b. Code State 


For each input bit to the encoder we have two output bits as shown in 
Figure 7. The number of states NS is. 


K-l 

NS = 2 = 32 


The rightmost five bits of the shift register constitute the state of the code. 
It is clear that the input bit will be the least significant bit for the next 
state. 


c. Quantization 

B 

We may quantize received samples into B bits. Thus we need L = 2 levels 
of quantization. For the binary code symbols x £ (0,1), with the following 
conversion. 
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0+1 


1 + -1 

we can assign integers 0 through L-l to tho quantitation levels, aa shown in 

Figure 8. 


OUTPUT (QUANTIZATION LEVELS) 



Figure 8 . Uniform L-Level Quantizer Where TH is Upper Threshold of Quantization 
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Table 2 shows some typical values for TH assuming the noise samples are 


normalized to have a unit variance. 

Table 2. Typical Values for TH 


L 

TH ■ - 1)A* 

2 

0 

4 

1 

8 

1.7 

16 

2.3 

32 

2.7 

64 

3.1 

128 

3.5 

256 

3.8 


A is threshold spacing. 


d. Metric Quantization 

Figure 3 shows integer code symbol metrics for L-level quantization. 


LY 


QUANTIZATION LEVELS 
2 L-2 L-1 


CODE SYMBOL 


L-1 

L-2 

L-3 


1 

0 

0 

1 

2 


L-2 

L-1 


Figure 9. Integer Code Symbol Metrics for L-Level Quantization 


The branch metrics can be computed from Table 3. 
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Table 3. Branch Metrics Computation for L-Level Quantization 


Branch Index 
IB 

Branch Code Symbols 
*1 *2 

Branch Metric 
LB 

1 

0 

0 

2(L-1) - LY X - LY 2 

2 

0 

1 

L-l - LY 1 + LY 2 

3 

1 

0 

L-l ♦ LY X - LY 2 

4 

1 

1 

LT 1 * “2 


The IB's are the branch Indexes associated with the various distinct branch 
symbols. For example for a branch with code symbols 10, IB * 3. 

e. Unquantized Metric 

If we do not use quantization for the received samples, then the code 
symbol metrics can be computed from Figure 10, and the branch metrics are 
computed from Table 4. 


CODE SYMBOL 



Figure 10. Code Symbol Metrics for Unquantized Samples 
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Table 4. Branch Metrics Computation for Unquantized Samples 


Branch Index 
IB 


Branch Code Symbols 
*1 *2 


1 

2 

3 

4 


0 

0 

1 

1 


0 

1 

0 

1 


Branch Metric 
LB 


r l ♦ *2 
r l - r 2 

~ r l ♦ r 2 

" r l - *2 


f . Vlterbi Algorithm 

The Vlterbi algorithm is essentially a maximum likelihood decoding pro- 
cedure for convolutional codes. The Viterbi algorithm recursively determines a 
path sequence with largest metric (most likely data sequence) through the trel- 
lis of the code from the initial state to the final state. The code trellis is 
basically an extended version of the code state diagram in time. The algorithm 
at any time retains the largest metric path to each state (node) and eliminates 
other smaller metric paths to the state from further consideration. The 
retained data sequence with the largest metric to each state Is called the sur- 
vived data sequence and the corresponding metric is called the survived metric. 
To implement the procedure, we must store the survived metric and the sur- 
vived data sequence for every state at any given time k. Then at time k ♦ 1, 
after computing the branch metrics, we can use the existing information which 
has been stored at time k, to compare and select the best metrics and the path 
sequences for each state and store them as the new survived metrics and the sur 
vived data sequences. We can proceed with this recursion, starting at time 
k - 0 at some Initial state, and ending the process at final time k ■ N. 
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Finally we ehooaa tha data sequence with largest metric among all atataa at time 
k ■ N as tha dacodad data sequence. 

In praetiea sinca N is usually large, tha amount of storage naadad to 
ratal n tha survivad data saquancas is larga. Tharafora if N is vary large, it 
is nacassary to truneata survivad saquancas to some langth a. However it can ba 
shown that there Is a high probability that all survived data sequences at tine 
k will have identical data bits very far back from the present tine k. This 
suggests that if tha algoritha stores enough of tha past data bits of each of 
tha 32 survivad sequences , than tha oldest bits on all stored data saquancas 
will ba identical. Our sinulation has shown that wa only need to store m » 36 
most recant bits of tha survivad data, sequences; in this case tha affect on per- 
formance is negligible. Than tha algorithsi outputs tha oldest bit in tha sur- 
vivad data sequence at state number 1, as tha hard dacodad bit. Another practi- 
cal consideration is metric overflows. Sinca the survivad nmtrics grow in time, 
to prevent a possible overflow, it is nacassary to renormalize the survivad 
metrics from time to time. As follows we explain the algorithm in detail. 


1. Definitions: 


V 1 J) 

LB (IB) 


LI 

V J > 


Survivor metric at state number J at time k. 

Array containing tha survivad metrics at time k, i.e., 
-k“ <M k < l >,M k <2> M k <32) > 

Branch metric corresponding to the branch coda symbols 

x, and x . indexed by IB 
12 

Array of branch metrics, i.e., LB • 

(LB(1) ,LB(2) ,LB(3) ,LB(4)) 

State number J; J ■ 1 , 2, . . . . , 32, at time k 
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MM 


k+1 


Array of the a aott raeont bits of tho survived data 
sequence, taralnatlng at state J at tlaa k 
36 x 32 aatrlx that stores all arrays 36 (J); J ■ 1,2, 
...» 32 at tlaa k 

3>*« x 32 a*trlx that storas all arrays 36 (J) ; 

J • 1,2 32 at tlaa k+1 


I& 

LI 


usr 


LT 


Array of all branch indaxas froa stata nuabar J to stata 

nuabar 2J - 1, for J • 1, 2, 16 

Array of all branch Indaxas froa stata nuabar 16 ♦ J to 

stata nuabar 2J - 1, for 3 • 1, 2 ,16 

Array of all branch Indaxas froa stata nuabar 16 ♦ J to 

stata nuabar 2J, for J ■ 1, 2, 16 

Array of all branch indaxas froa stata nuabar 3 to stata 
nuabar 2J, for J - 1, 2, 16 


Nota that sinca tha first and tha last stages of tha shift raglstar ara 
connactad to tha MOD 2 addars, than 


]&' - Lfi and LT ■ LI 


2. Stata Transitions 

For a givan J; 1 < J < 16, tha stata transitions ara shown in 
Figure 11, where IB ^ IB'. All state transitions associated with the branch 
Indexes are shown in Table 5. Note that when we use branch indexing we need 
only to compute 4 metrics each time rather than computing all metric branches. 
Therefore branch metric computation is almost independent of number of states. 
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Figure 11. State Transitions Between Pair of States 


l 
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Table 5. State Transitions With Associated Branch Indexes 


8 k (J) 

S k+1 <2J-l> 

XB 

1 

1 

1 

2 

3 

4 

3 

5 

3 

4 

7 

2 

5 

9 

1 

6 

11 

4 

7 

13 

3 

S 

15 

2 

9 

17 

3 

10 

19 

2 

11 

21 

1 

12 

23 

4 

13 

25 

3 

14 

27 

2 

15 

29 

1 

16 

31 

4 


i 

* 

i 

> 
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Tabic 5. State Transitions With Associated Branch Indexes (Continued) 


S k <16+J) 

S k+ i<2J-l> 

IB 

17 

1 

4 

18 

3 

1 

19 

5 

2 

20 

7 

3 

21 

9 

4 

22 

11 

1 

23 

13 

2 

24 

15 

3 

23 

17 

2 

26 

19 

3 

27 

21 

4 

28 

23 

1 

29 

25 

2 

30 

27 

3 

31 

29 

4 

32 

31 

1 
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Tablt 5. State Tranaitlona With Associated Branch Indexaa (Continued) 


i k (J) 

3»fl<2J> 

IB 

1 

2 

4 

2 

4 

1 

3 

6 

2 

a 

8 

3 

5 

10 

4 

6 

12 

1 

7 

14 

2 

8 

16 

3 

9 

18 

2 

10 

20 

3 

11 

22 

4 

12 

24 

1 

13 

26 

2 

14 

28 

3 

15 

30 

4 

16 

32 
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Table 5. State Transitions With Associated Branch Indexes (Continued) 


S k (16+J) 

S k+1 (2J) 

IB 

17 

2 

1 

18 

4 

4 

19 

6 

3 

20 

8 

2 

21 

10 

1 

22 

12 

4 

23 

14 

3 

24 

16 

2 

25 

18 

3 

26 

20 

2 

27 

22 

1 

28 

24 

4 

29 

26 

3 

30 

28 

2 

31 

30 

1 

32 

32 

4 
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3. Algorithm 


Step 1. Stores*: 

Store Branch index arrays, 

LO - (1, 4, 3, 2, 1, 4. 3, 2, 3, 2, 1, 4, 3, 2, 1, 4) 
LI - (4, 1. 2, 3, 4, 1, 2, 3, 2, 3, 4, 1, 2, 3, 4, 1) 
Store the survived data sequences in buffers, 

"SlU ' <! Vh, 36 (1) ’ “k+l,36 <20 \+l,36 <32)) 

Store branch metric array, 

LB - ( LB ( 1 ) , LB ( 2 ) , LB(3), LB(4>) 


< \.36 <1> - W 2> W 32 ” 


Step 2. Initialization: 

Buffers, 

MO * (0,0, .... 0) 

MN(1, 2J-1> * 0 

I J * 1, 2, 16 

MN(1, 2J) - 1 
Metrics, 

M 0 (13) * 0 

M 0 <I) « -2000 ; I - 1, 2, 3, ...» 32, I # 13 

Stating time, k « 0 


Step 3 Recursion: 
k - k + 1 

branch metric computation, 

compute branch metrics from Table 4 and store in LB 
I - 0 
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Step 4 I - I + 1 

Survived metrics, 

\ (21-1) - Max | fM k _ 1 (D + LB(LO(I)) J, 

[M kl (16+I) + LB (LI (I)) ] | 

M k (21) - Max | [M k _ 1 (I) + LB(L1(I))J, 

[M kl ( 16+I) + LB(LO(I)) ] ^ 

Select the survived sequences from M0 R ^ , add selected bit "0" 

or "1" to it, store in scratch buffer MN R , 

If M k (2l-1) = M k _ 1 (I) + LB(L0(I)) 

then MN k (J+l, 21-1) = M0 R _ 1 (J, I) ; J - 1 , 2 , . . . , 35 

If M k (2I-l) * M k _ 1 (16+I) + LB (L0(I) ) 

then MN R (J+l, 21-1) = M0 R _ 1 (J, 16+1) ; J = 1, 2, .... 35 

If M k (2l) = 14^(1) + LB(LIU)) 

then MN R (J+l, 21) - MO kl (J, I) ; J * 1, 2, .... 35 

If M k (2I) = M^I) + LB(LO(I)) 

then MN R (J+l, 21) = M0 R _ 1 (J, 16+1) ; J = 1, 2, . . . , 35 

If I < 16 go to step 4 

restore MN, into MO, 
k k 

M0 R ( I , J ) = MN k (I, J) ; I = 1, 2, 36 

J = 1, 2, .... 32 

The decoded bit at time k is 

MO R (36, 1) 
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Normalization of survived metrics, 


Find the largest M^J); J - 1. 2. .... 32 end call it MXM, then 
M^J) - M^J) - MXM j J - 1, 2, .... 32 
call node synchronisation algorithm and store the result in NOD. 

If K < N go to step 3 

Otherwise, choose the largest metric aawng all the survived metrics at 
k ■ N. The corresponding data sequence in the buffer MO is considered 
as the final decoded sequence. 
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g. Bit Error Rate Performance 


From our analysis using Table 1 for weight distribution of code and then 

£ 

using union bound we can approximate [3] the bit error rate P, as 

b 



where 


— is bit signal to noise ratio and 

N o 


erfc(z) 


j_r .t= 
vrj 


dt 


Figure 12 gives the theoretical performance of this algorithm and the 
simulation results. The .25 dB loss in simulation with respect to theoretical 
result is due to infinite bit quantization that we assumed in theoretical result, 
and it is due to path memory truncation and metric quantization in the software 
algorithm. 


theoretical upper bound is [2] 



where a, are coefficients in tr 'nsfer function bound [2]. The first few 
d 

are given in Table 1. 
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IX. CONCLUSION 


The algorithms and subroutines developed in this study are suggestions for 
constructing an actual data processing program and may require some further 
optimization. 

Simulations of the carrier tracking loop have indicated that the least- 
squares procedure for finding w and w may not converge if the starting parameters 
are not close enough to the actual w and w, if the noise is too large, or if the 
phase fluctuations are so large that the 30 second coherence span is too long. 
Normally if the carrier extraction is successful, the subcarrier extraction, bit 
and frame synchronization proceed smoothly. Simulations of the decoding algo- 
rithm are in good agreement with the theoretical performance. Detailed analyses 
of the statistical performance of the subcarrier demodulator has not been carried 
out, but the techniques used for signal processing should not produce any signi- 
ficant degradation from the theoretical values. For this reason, the major un- 
certainty in the decoding is associated with the carrier extraction and the 
processes that could cause unmodeled phase modulation of the carrier wave. 
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APPENDIX 


SUBROUTINE SETUP 

C***************** SET UP FOR VITERBI DECODER **************** 

C M IS VECTOR OF STATE METRICS 

C LO IS INDEX VECTOR OF ALL STATE TRANSITIONS SUCH AS J->2J-1 
C » J*l» 2* * 1 A. 

C LI IS INDEX VECTOR OF ALL STATE TRANSITIONS SUCH AS 
C 16*J->2J-1j J»l,2, , 16. 

C NOTE : INDEX VECTOR OF ALL STATE TRANSITIONS SUCH AS J->2J 
C IS LI. 

C INDEX VECTOR OF ALL STATE TRANSITIONS SUCH AS 

C 16*J->2J IS LO. 

C MO IS 36*32 MATRIX WHICH STORES ALL BEST PATHS TO 32 STATES, 

C IT ONLY KEEPS PATHS WITH LENGTH 36. THE 37TH BIT CAN BE 
C REGARDED AS HARD DECODED BIT. 

C MN IS 36*32 SCRATCH MATRIX. 
C*************************************************************C 
COMMON/LA l/LO< 16), LI < 16), M<32), MO<36, 32), MN(36. 32) 

DATA LO/1, 4, 3, 2, 1, 4, 3, 2, 3, 2, 1, 4, 3, 2, 1, 4/ 

DATA Ll/4, 1, 2, 3, 4, 1, 2, 3, 2, 3, 4, 1. 2, 3, 4, 1/ 

C INITIALIZATION 
DO 9 J-1,32 
DO 9 1-1,36 

9 MO < I , J ) -O 

DO 10 1-1, 16 

11- I+I-l 

12 - 1*1 

MN( 1 , ID-0 

10 MN< 1, 12) — 1 

C ASSUME ENCODER IS IN STATE « 13 
DO 20 I =>1,32 
20 M<I)— 2000 

M( 13) -O 
RETURN 
END 


SUBROUTINE DECODE (KTMI, LY1, LY2, L, KTMO, KOUT) 

THIS SUBROUTINE IS VITERBI DECODING ALGORITHM 
KTMI IS INPUT TIME 
KTMO IS OUTPUT TIME 
LY1 AND LY2 ARE QUANTIZER OUTPUTS 
L IS NUMBER OF QUANTIZATION LEVELS 
KOUT IS DECODE BIT, IF IT IS NOT EQUAL TO 9 
MAXN IS THE STATE NUMBER WITH LARGEST SURVIVED METRIC 
DIMENSION LB ( 4 ) , MS ( 32 ) 

COMMON/LA 1 /L0( 16 ) , LI ( 16 ) , M( 32 ) , M0< 36, 32),MN<36, 32) 

COMPUTATION OF BRANCH METRICS 

i 

i 


A-l 


nooooooooooooo 


c 


CALL BRANCH <LY1,LY2, LB,L> 


C SELECTING THE SURVIVED METRICS AND SURVIVED DATA SEQUENCES 
DO 10 K-l, 16 
I J-K+K— 1 
XL-K+K 
XK-K+16 

MTO-M(K>+LB<LO<K>> 

MT1«M(XK>+LB(L1<KM 
MT2-M<K)+LB<L1<K> > 

MT3— M( IK)+LB(LO(K) > 

IF(MTO. OE. MT1 ) QO TO 20 
MS ( I J ) — MT 1 
DO 30 1-1,35 
Il-I+l 

30 MN<I1, XJ)-MO<I, IK) 

00 TO 40 
20 MS< I J)— MTO 
DO 30 1-1,33 
Il-I+l 

30 MN ( X 1 , I J ) -MO ( I , K ) 

40 CONTINUE 

IF(MT2. OE. MT3) 00 TO 60 
MS( ID-MT3 
DO 70 1-1,33 
Il-I+l 

70 MN(I1, IL)-MO(I, IK) 

00 TO 10 
60 MS( IL>— MT2 
DO SO 1-1,33 
I l-I+l 

80 MN(I1, IL)-MO(I, K) 

10 CONTINUE 
C 

C NORMALIZE M 
C 

DO 130 1-1,32 
150 M < I ) —MS ( I ) 

CALL MAX <MXM, MAXN) 

DO 90 1-1,32 
90 M< I )— M< I )— MXM 
C 

DO 100 J-l, 32 
DO 100 1-1,36 
100 M0< I, J)-MN< I , J) 

C OUTPUT TIME AND DECODED BIT 
KOUT-9 

IF(KTMI. OE. 36) K0UT*M0<36, 1 ) 

IF <KTMI. OE. 36) KTMO-KTMI-33 
RETURN 
END 
C 
C 
C 

SUBROUTINE BRANCH(LY1, LY2, LB, L) 

C THIS SUBROUTINE COMPUTES THE BRANCH METRICS 


A- 2 


non oooo non oo ooo 


C LY1 AND LY2 ARE QUANTIZER OUTPUTS 
C L 18 NUMBER OF QUANTIZATION LEVEL8 
C ARRAY LB STORES 4 BRANCH METRICS 
DIMENSION LB< 1 ) 

LB( 1 >-L+L-2-LYl-LY2 

LB ( 2 ) -L- 1 -LY 1 +L Y2 

LB<3)-L-H-LY1-LY2 

LB < 4 ) -L Y 1 H.Y2 

RETURN 

END 


SUBROUTINE MAX (MXM. MAXN) 

MXM IS LARGEST SURVIVED METRIC 

MAXN IS THE STATE NUMBER WITH LAROEST SURVIVED METRIC 
C0MM0N/LA1/L0(16).L1<16),M<32). MQ<36. 32) . MN(36. 32) 
MXM-M(l) 

MAXN-1 

DO 10 1-2.32 
IF(M(I). GT. MXM) 00 TO 20 
00 TO 10 
20 MXM-M(I) 

MAXN— I 

10 CONTINUE 
RETURN 
END 


SUBROUTINE QUANT <R. L. LY. TH) 

THIS SUBROUTINE IS L LEVEL QUNTIZER 
TH IS UPPER LIMIT OF QUANTIZATION 
R IS RECEIVED SAMPLE AT THE INPUT OF QUANTIZER 
LY IS OUTPUT QUANTIZATION LEVEL 
SPACE- < TH+TH ) / ( L-2 ) 

IF(R. OT. TH) LY-0 
IF(R. LE. TH) LY-L-1 
IF(ABS(R). LE. TH) 00 TO 10 
00 TO 20 

10 LY-(TH-R) /SPACE* 1 
20 CONTINUE 
RETURN 
END 
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